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Abstract 

Combining an electronic theory with molecular dynamics simulations we 
present results for the ultrafast structural changes in small clusters. We de- 
termine the time scale for the change from the linear to a triangular structure 
after the photodetachment process Ag^ — > Ag 3 . We show that the time- 
dependent change of the ionization potential reflects in detail the internal 
degrees of freedom, in particular coherent and incoherent motion, and that 
it is sensitive to the initial temperature. We compare with experiment and 
point out the general significance of our results. 
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The excitation of a cluster by a laser pulse induces time-dependent changes of its elec- 
tronic and atomic structure. These changes involve bond formation and bond breaking. The 
understanding of the relaxation mechanisms induced by the excitation is of general interest. 
In particular, it is of fundamental importance to study how the system approaches equilib- 
rium in order to determine how the time-scales of the relaxation processes can be controlled 
by varying the experimental conditions. 

Recently, a pump&probe experiment has been performed on mass selected Ag3 clusters, 
which serves as an example for the investigation of the structural relaxation times [l]. The 
initially negatively charged clusters were neutralized through photodetachment by the pump 
pulse and after a delay time At ionized by the probe pulse in order to be detected. Due 
to the remarkable differences in the equilibrium geometries of the ground states of Ag3 , 
which is linear 0, and Ag 3 , which consists of an obtuse isosceles triangle H, the ultrashort 
photodetachment process puts the neutralized trimer in an extreme nonequilibrium situa- 
tion. As a consequence of that, a structural relaxation process occurs. The experimental 
signal, consisting in the yield of Ag^ was measured as a function of At and the frequency 
of the probe laser pulse. For a frequency slightly above the ionization potential (IP) of Ag3 
a sharp rise of the signal is observed at At ~ 750fs. After a maximum is reached, there is a 
saturation of the signal, which then remains constant for at least lOOps, which is the longest 
time delay used in the experiment. New features appear for higher frequencies. Again the 
signal increases sharply, but after reaching a maximum it decreases to a constant value. 

A preliminary interpretation of these results uses the Franck-Condon-principle Q. The 
first laser pulse creates a neutral linear silver trimer which bends and comes to a turning 
point near the equilateral equilibrium geometry of the positive ion ||. After rebounding, 
the neutral trimer starts pseudo-rotating through its three equivalent obtuse isosceles equi- 
librium geometries. This would explain the saturation behaviour of the signal. However, 
this would mean that the pseudo-rotations have an extremely long mean life, which seems 
improbable. Furthermore, this model does not explain why the signal changes as a function 
of the frequency of the laser pulse. 
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In this paper we perform a theoretical analysis of the physics underlying the ultrafast 
dynamics of Ag 3 clusters produced by photodetachment. In particular, we analyze the time 
evolution of the ionization potential and the dependence of the dynamics on the initial 
temperature of the clusters. We show that the experimental results can be explained using 
a physical picture which can be generally applied to other ultrashort-time processes. In our 
calculations, we combine molecular dynamics (MD) simulations in the Born-Oppenheimer 
approximation with a microscopic theory to describe the time-dependent electronic structure 
of the clusters. 

In order to determine the potential energy surface (PES) needed for the MD simulations, 
we start from a Hamiltonian of the form H = HTB + ^/2J2i^j 'K 1 "*?); where the tight-binding 
part Htb is given by 

HtB = £ ia C taa C iacr + X! ^aj'/3 C ia C i/3- (1) 

i,a,cr j^j,<7 
a,f3 

Here, the operator cf aa (cj Q<T ) creates (annihilates) an electron with spin a at the site % 
and orbital a [a = 5s,5p x ,5p y ,5p z ). e ia stands for the on-site energy, and V ia jj3 for the 
hopping matrix elements. For simplicity, and since the 5s electrons are expected to be 
rather delocalized, we neglect the intraatomic Coulomb matrix elements. <p{^ij) refers to 
the repulsive potential between the atomic cores i and j. For the distance dependence of 
the hopping elements and the repulsive potential we use the functional form proposed in 
Ref. U|. By diagonalizing H T b, (taking into account the angular dependence of the hopping 
elements 0), and summing over the occupied states, we calculate as a function of the atomic 
coordinates the attractive parts of the electronic ground-state energies E~ ttr and E® ttr of Ag 3 
and Agg, respectively. Then, by adding the repulsive part of H we obtain the PES, which 
we need to perform the MD simulations. In order to determine the forces acting on the 
atoms we make use of the Hellman-Feynman theorem. Thus, the a-component of the force 
acting on atom % F ia = —dE/dr ia is given by 

^ = -E<*lf^l*>-iE^- (2) 

occ. 
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Here the | k)'s are the eigenstates of Htb- The parameters of H are determined in the fol- 
lowing way. The on-site energies were obtained from atomic data || . We fit the parameters 
V a f3 and the potential (j>{^ij) in order to reproduce the equilibrium bond lengths of the Sil- 
ver dimers Ag2 , Ag2 and Ag^ obtained by effective core potential-configuration interaction 
calculations [@,§J. We assumed the hopping elements V a p to fulfill Harrison's relations ||. 
The best fit was obtained by the following parameters: V sp = 0.954eV, r c = 4.33A, n c = 2, 
m = 5.965 and A = 0.605eV, where r c and n c refer to the cutoff radius and exponent, 
whereas m and A stand for the exponent and strength of the repulsive potential <ft(rij) @. 
Using these parameters we have calculated the vibrational frequencies of the dimers, which 
compare reasonably well with the experimental values |IJ and quantum-chemical calcula- 



tions @3- Then, we determined the equilibrium geometries of the ground states of the silver 
trimers Ag 3 , Ag 3 and Ag^ which again yielded excellent agreement with the all- valence elec- 
tron calculations of Refs. and 0. We have also determined a linear equilibrium geometry 
for the excited state 2 S+ of Ag 3 , which can be reached by photodetachment of Ag 3 ||. 
However, in our calculations this state lies 1.2eV higher than the ground state (in relatively 
good agreement with photodetachment experiments ||) and therefore cannot be excited by 
the laser frequencies used in Ref. [|I). Therefore we consider in the following only the dy- 
namics of the ground state of Ag 3 . Note, our hopping parameters are comparable to those 
of silver bulk [^]. This gives another justification for the neglect of Coulomb interactions. 

The MD simulations are performed applying the Verlet algorithm in its velocity form. 
We used a time step of At = 0.05/s. This ensures an energy conservation up to 10 _6 eV 
after 10 5 time steps. The equilibrium structures were obtained by performing simulated 
annealing. Starting with the equilibrium geometry of Ag^, we generate an ensemble of 
approximately 1000 clusters characterized by the ensemble temperature T, defined as the 
time average of the kinetic energy for a long trajectory (~ 10 6 time steps) [ [TT|j . 

In Fig. 1(a) we show the time dependence of the fraction of Ag 3 clusters p{hv,t) with 
ionization potentials IP(t) < hv. We scale hv with IP , which is the minimal ionization 
potential of Ag 3 in the equilateral equilibrium geometry of Agg". The quantity p(hi>,t) 



can be interpreted as the time- dependent probability for ionization with a laser pulse of 
frequency hv (i.e., half of this energy, if the ionization is achieved in a two photon process as 
in Ref. [J]). p(hu,t) is proportional to the signal detected in the experiment. Note, slightly 
above the ionization threshold (Jnv = 1.02IPo), clusters can be ionized only after a delay time 
to(hi') ~ 750fs, at which p(hv,t) begins to increase. p(hu,t) displays the features observed 
in the experiment, namely first a sharp increase, then a maximum and finally a decrease to 
a smaller value, which remains constant. There is, however, a narrow range of frequencies 
hi>, where the enhancement of the signal after the sharp increase is largest. For higher laser 
frequencies the maximum becomes broadened and the enhancement factor is smaller. 

These results can be understood physically as follows. Upon photodetachment of a 
binding electron vibrational excitations occur, in particular, those of the central atom along 
the chain direction. Due to the shape of the PES the motion of the central atom dominates 
the ultrashort time response over the first few hundred fs. Then, the slower thermally 
activated bending motion comes into play and yields triangular bonded Ag 3 . The resultant 
bond formation is exothermic. The excess energy can in turn cause bond breaking or in case 
of uniform energy distribution also a regular vibrational mode like pseudo-rotations. 

In order to demonstrate the temperature dependence of the bond formation and bond 
breaking processes, we present in Fig. 1(b) results for the time dependent probability 
Vtriangif) of finding triangular clusters for different initial temperatures. Ptriang{t) is de- 
fined as the fraction of clusters having all angles between 45 deg and 90 deg. Comparison 
of the behaviour of ptriangif) for T = 317K with the curves of Fig. 1(a) clearly shows that 
only these particular cluster geometries contribute to the signal observed in the experiment 
Jl[. Note that the onset of the abrupt increase of ptriangif) shows a strong temperature 
dependence. For increasing initial temperature the clusters bend faster. Notice also the 
remarkable fact that at low temperature the maximum in pt r iang{t) disappears. 

Fig. 2 shows contour plots for the time development of the distribution of the IP's of 
neutral silver trimers after photodetachment at t = 0. The most prominent feature of these 
pictures is the coherent oscillation of the IP of all clusters, which continues for about 700fs 
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in the case of clusters starting at T = 29K and which is destroyed much more quickly at the 
higher temperatures. These coherent oscillations of the IP are due to the internal vibrations 
of the linear chain before and during bending. High values of the IP correspond to a rather 
symmetric geometry, while the lower values stand for the turning points of the oscillation. 
The second feature to be noted are the areas in the plot with IP's below approximately 
1.07IPo, that become populated only after time delays of 1500, 700 and 300 fs respectively 
and which represent those clusters that have bent to form geometries close to the equilibrium 
geometry of Ag 3 . For T = 1048K we observed fragmentation of 31% of the trimers (after 
20ps), and in Fig. 2(c) we show only the distribution corresponding to the non-fragmented 
clusters. 

With the help of Fig. 2 it is also possible to explain the maximum of the IP as a result 
of the coherent vibrational motion of the clusters: while a large fraction of the clusters 
bend collectively, they still continue their vibrational motion, which results in a pronounced 
minimum of the IP in the same way as for the unbent trimers. This effect accounts for the 
dark areas in Fig. 2(b) and 2(c) for IP(i)/IP below 1.07 approximately and is quantified 
by integration over energies in the calculation of p(hu,t) [Fig. 1(a)]. Since the clusters 
responsible for this effect have gained an average energy of 900K while descending the PES 
towards the equilibrium structure of Ag 3 , their coherent motion is quickly destroyed and 
this explains why no further maxima occur. This is supported by the fact that the cluster 
ensembles starting out with temperatures of 29K and 104K do not show this maximum. 
Here, the coherent motion has already disappeared before the clusters bend toward the 
equilibrium geometry of Ag 3 . 

Fig. 2 reveals that for longer times (t > 2000fs) the clusters approach a time-independent 
distribution over the IP's or equivalently over the Ag 3 geometries. These asymptotic distri- 
butions are plotted in Fig. 3(a). 

In Fig. 3(b) we have plotted the quantity p(hu,t) as in Fig. 1(a), but for an initial 
temperature of T = 29K and for a higher energy hv = 1.17IPo- A new feature appears 
which consists in sharp oscillations of the probability between and 1 within the time range 
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[0 < t < 700fs] . These oscillations result from the dominant vibrational mode of the trimers 
for the first 700fs approximately. This interesting effect was not observed in the experiment 
|Ij], because the laser energies used were not high enough. In order to demonstrate that the 
excited state 2 S+ does not contribute to the signal observed in the experiment of Ref. we 
have performed MD simulations on the corresponding PES and confirmed its one dimensional 
dynamics. Its signal p(hu, t) consist in oscillations with no sharp increase like that observed 
in Ref. |lj]. However, it is important to point out that this structure-less signal could interfere 
with the bending or stretching dynamics of the ground state if the state 2 £+ were excited 
by the pump pulse. 

The last important point to be investigated is the role of pseudo-rotations. In order to 
study these vibrational excitations we determine the autocorrelation function 

G{r) = (3 ) 

z2k Qx\tk)Qx{tk) 



where q x is one of the normal coordinates of the triangular Ag 3 molecule [12] . If the pseudo- 
rotations had a long life time, G(r) would show an oscillatory behaviour. However, as 
shown in the inset of Fig. 3(c), the correlations are strongly damped. This becomes clear by 
analyzing the Fourier transform of G(t). The life time of the pseudo-rotations, determined 
from the width of the peak at half maximum, is only slightly lower than the duration of 
one cycle. This shows that pseudo-rotations do not play an important role in the dynamics 
of the trimers. The high kinetic energy of the atoms seems to prevent the occurrence of a 
regular mode like pseudo-rotations. The saturation of the signal is rather a statistical effect 
induced by the temperature. 

In summary, by employing an electronic theory and MD simulations we have analyzed 
the femtosecond dynamics of Ag 3 upon photodetachment. Our results show that an uncom- 
plicated theory which can easily be extended to deal with larger systems is able to account 
for all experimental observations. Most importantly, we find that the time scales for bond 
formation and breaking is affected by temperature. Our results might also be of general 
interest regarding laser control of chemical reactions. Fig. 3(a) suggests that temperature 



can be used to design special distributions for the IP or for the cluster geometries present 
in the initial ensemble. Fig. 3(b) immediately suggests new experiments with a probe laser 
pulse energy of hu = 1.17IPo- Thus a direct measurement of the dominant vibration of the 
trimers seems possible. 

This work was supported by the Deutsche Forschungsgemeinschaft. 
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FIGURES 

FIG. 1. (a) Time dependence of the fraction of clusters p(hu,t) with IP(i) smaller than 
his. The sharp increase and the overall time dependence of p(hv, t) for increasing his should be 
compared with the experimental results of Ref. 1. The initial temperature was T = 317K. (b) 
Time dependence of the fraction of clusters ptriang(t) having a triangular structures for different 
initial temperatures. 

FIG. 2. Distribution of photodetached clusters as a function of time and IP at different 
temperatures. Black regions indicate large number of clusters, whereas no clusters are present in 
the white parts. Note the progressive destruction of the coherent motion after photodetachment 
for increasing initial temperatures. 

FIG. 3. (a) Average distribution of neutral silver clusters over IP's for 2ps < At < 20ps 
after photodetachment. (b) Calculated p(hu, t) for a large value of his. Note that also the internal 
vibration of the chain before and during bending can be detected, (c) Spectral analysis of the 
pseudo-rotations. Note that the life-time is smaller than one pseudo-rotation period. 
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